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. Abstract 


The quantitative interpretation of infrared sunset spectra in terms 
of mixing ratio profiles of atmospheric constituents requires a 
detailed consideration of the airmass and pressure distribution 
along the sunset ray. A computer program has been developed for 
determining the airmass and pressure distribution along this ray. 
Since infrared absorptions are pressure-dependent, a technique has 
been developed for treating this pressure dependence in analyzing 
absorption over such paths. The technique uses an approximation 
which maintains computational accuracy for such paths while reducing 
the computer time required for the calculations. 
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ABSTRACT 


The quantitative interpretation of infrared sunset spectra in 
terms of mixing ratio profiles of atmospheric constituents requires a 
detailed consideration of the airmass and pressure distribution along the 
sunset ray. A computer program has been developed for determining the 
airmass and pressure distribution along this ray. Since infrared absorp- 
tions are pressure-dependent, a technique has been developed for treating 
this pressure dependence in analyzing absorption over such paths. The 
technique uses an approximation which maintains computational accuracy 
for such paths while reducing the computer time required for the calcula- 
tions. 

The interpretation of such spectra is further complicated if the 
concentration of the constituent of interest varies with solar zenith angle 
(e.g.NO). Calculations were performed using time-varying profiles for 
NO and NO.. These calculations show that, if this time variability is not 
taken into account, the inferred profiles are greatly in error in the case 


of NO and marginally in error for NO,- Kecession For 
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I, Introduction 


A powerful infrared technique for the detection of trace consti- 
tuents is to observe atmospheric transmission spectra from an airborne 
platform using the sun as a source. The sensitivity of the technique 
is enhanced by obtaining the spectra at solar zenith angles >90° 
(sunrise or sunset). In fact, within the last decade, several trace 
constituents have been measured in the lower stratosphere by this 


technique: i.e., HNO,, NO,, NO, HCl, HF, CF,Cl,. cFCl, and 
ccl 


: The quantitative derivation of mixing ratio altitude profiles 
from sunset solar spectra requires a detailed consideration of the 
geometry of the optical path through the atmosphere, including the 
effect of atmospheric refraction. 

The problem is further complicated by the fact that infrared 
absorptions depend on the pressure and temperature environment of 
the molecular species responsible for the absorption. In addition, 


when obtaining profile data by analyzing spectra obtained at solar 


zenith angles greater than 90°, it is generally assumed that the 


altitude distribution of the constituent of interest does not change with 
the time and geographic location of the optical path during the time 
required to obtain the spectra. However, these assumptions are not 
valid for photochemically-active species (e.g., NO). The objective 

of this program has been to investigate the importance of these effects 
in the analysis of sunset solar spectra. Computer programs have 
been developed that accurately determine the airmass traversed by 
solar radiation in each portion of the optical path. In addition, a 
technique has been developed for including the pressure and temperature 
dependences of the absorption that does not use an excessive amount 
of computer time. Finally, these programs have been used, along 


with time-dependent profiles for NO and NO,, to investigate the errors 


1 


in the profiles which are introduced when time dependence is not a 


part of the spectral analysis. The time-dependent profiles were supplied 
by Dr. Megill of Utah State University under sub-contract. Detailed 

discussions of the various aspects of the problem, which are summar- 
ized here, have been presented in reports issued during the stegrens: 


and are listed below. 


Il. Airmass Distributions 


The geometry for the problem of low sunset data is shown 
in Figure 1. The major point of interest is that the concentration 
along the path of a photochemically-active constituent is quite different 
from the sun to the tangent point (minimum altitude) than it is from 
the platform to the tangent point, since the solar flux is different on 
the two sides. The airmass distribution along a sunset path is derived 


using a ray-tracing program described in detail in previous reports 


[Snider and Goldman, 1975; Snider, 1975]. Calculations were performed 


using, as input, a standard midlatitude summer atmosphere [McClatchey 
et al., 1972] and 5ym as the wavelength of interest. Typical results 

of such airmass and accumulated airmass calculations are shown in 
Figures 2 through 5, for an observer at 30 km and apparent zenith 
angles of 92° and 94°, Calculations for other altitudes and zenith 
angles were also performed and these figures and tables are included 
in Goldman et al. [1977]. Table I is an example which shows the 
results in tabular form of such a calculation for 30 km and 94°. The 


tabulated effective pressure, P .., and effective temperature, T 


eff eff’ 
are average values over the 1 km shells used in the program as 


2! he cloak 
Peft — a cig ee <i (1) 
dm dm 


where dm is the differential airmass. 
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The integrands for these experiments are given in the tables as 
DEL(P-DM) and DEL(T-DM) and the airmass increments are listed 

as DEL(DM). The values are given for the two sides of the path: 
from the tangent point to the observer and from the tangent point 
toward the sun. The angle given is the zenith angle at the intersection 
of the path and the boundary of the shell. Finally, the geometrical 
path and airmass are given as ACCUM DEL(PATH) and ACCUM DEL(DM). 


Ill, Interpretation of Infrared Sunset Solar Spectra 


A. Weak line approximation 


Denote the observed integrated absorption at different apparent 


zenith angles z Zz Been le5y Ay Ay, -.., for a given observer 


altitude, h,- ona a loss of generality, assume that h, is an 
integer km altitude. For simplicity of the discussion, the weak 

line approximation is assumed, so that the integrated absorption, A, 
is given by hie) = Steui’  Triahacdla oe} U luteculefon), 
where S is the integrated line intensity and U is the absorber amount. 
The airmass distribution tables [Table I here and those in Goldman 

et al., 1977] give atmospheric amounts U as ina M is the F 
tabulated airmass, and Uj= 1 airmass = 2.153 x 10° molecule/cm . 
Trace constituent amounts are described by BU where 8 is the gas volume 
mixing ratio. 

The tables are used for a layer-by-layer derivation of a minor 
constituent altitude profile in the following way. The atmosphere is 
divided into even 1 km shells, although Figure 1 shows 2 km shells 
for graphic simplicity. The 1 km shells listed in the tables can be 


accumulated over two or more km as required. Starting with the 


first path and Zy > 90° (required to observe the weak absorption of 


the trace constituent), a constant mixing ratio, B,, can be derived 


1’ 


See 


9 


PO a aaa Ere 


(1) 


for the atmosphere above Fain of the path according to 


By = SUM, (2) 


where M, is the total airmass along path 1. This B, will be con- 
sidered as the mixing ratio for layer 1, between h, and h,, as well 
as for the residual atmosphere above h, in the sun direction. It 
should be noted that H ain can obviously fall between the even km 
boundaries, so that h, will be the same altitude as A nin or the next 
nearest even km below it. Subsequent inner layer boundaries are 
defined by the minimum altitudes of the corresponding zenith angles, 
Z » 80 that a layer can include more than 1 shell with increasing 
zenith angle. For example, with the observer at 30 km, the layer 
defined by z = 92° and z = 93° includes five 1 km shells. This is 
shown as Layer 3 on Figure 1. 


A mixing altitude profile is then derived by the successive 


equations: 
n-l 
“ (n) (n) 
of alee ds ee ee, ee (3) 
6M sum ,™) 
o n,ntl 
Here, BL is the constant mixing ratio for the layer defined by H sig 


(n-1) (n) 
1 min ~ k+l 
portion of the airmass between altitudes hy, and Ley for path n. The 


to ho » the even km boundary below H 
upper index n denotes the different paths for a given observer altitude. 
The M,. el (n) values are those tabulated in the airmass distribution 
tables every 1 km, accumulated over the proper layers. In Eq. (3) 
the summation in the numerator represents the correction to AL due to 


the absorptions in the outer layers, as each one is crossed twice. A 


10 


smaller correction is a™ = psu M_™) » where M,. is the residual 


airmass outside hy in the sun direction, The inner layer n is crossed 
once, but its airmass is twice that from H (n) - 
min n-l 


The corrections to A, are by no means negligible. However, 


a commonly used approximation is to neglect these corrections in 
(n) 

ntl 

total airmass along path n, and assign the resulting mixing ratio, 


the numerator and replace M, in the denominator by M.: the 
’ 


A, 
B n iy SU M (4) 
on 
to H , or a small altitude range above it. This approximation is 


motivated by the fact that a significant portion of the total airmass 
along the path is indeed concentrated near Hain’ It can be a good 


approximation only if the effect of neglected corrections in the numera- 


(n) 


tor are compensated for by the ratio M/ M, nil so that 
-l M 
(n) (n) n ba: 
(a, z Cal 1 BM. k+l - 4 ) (n) 4, : 


n, ml 


Examination of the airmass distribution tables shows that for 
many practical cases, this approximation can lead to errors from 5% 
to a factor of 5. 

Another commonly used approximation is to consider the absorp- 


tion difference, AL - AL 1’ and the corresponding total airmass difference, 


M_ - M 
n n- 


1° 8° that B. is approximated by 


(n) 


and assigned to H nin as in the previous approximation. This 


approximation can be quite good if both A, > AL and M, > M, 


-1 
and M, is not very large 


oe kg 


but in many practical cases AL ws AL 1 


compared to M unless z 294°, A combination of these conditions 


i) 
n-l 
can result in a very large error. 


B. General Case 


The layer-by-layer method described above can be extended 
to a more general line-by-line layer-by-layer calculation over a 
spectral interval of non-weak overlapping lines. This can be done by 
an extension of the method presented by Goldman et al. [1973a; 
1973b], to be referred to hereafter as I and II. The extension 
presented here is for both small and large zenith angles. 

The mixing ratio is derived by fitting the theoretical line-by-line 
integrated absorption to the measured value for successive atmospheric 
layers, starting from the uppermost layer. For each ray path, 
except for the first one, the equations used involve only two successive 
layers, but with an adjustment for the changing direction of the path. 
For high sun, as occurs for example during. balloon ascent, the adjust- 
ment is based upon the secant of the zenith angles of the different rays. 
For low sun, as occurs in sunset observations from balloon float 
altitude, the adjustment is based upon the actual airmass distribution 
along each path. In both cases, the contribution from all layers along 


the path is taken into account. 


1. High Sun 


The geometry for high sun data is shown in Figure 6, assuming 
that atmospheric curvature and refraction are negligible. The balloon 


altitude H Pp! corresponding to specific measurements, are 


prey» 


shown on one radial line, and the zenith angles Gis +s. O. increase 


(in this case) with altitude. At each altitude, 2 a ray n intersects 
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Fig. 6. 
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Solar ray paths for absorption measurements made for 
high sun at different altitudes during balloon ascent. 


13 


RAY | 


the radial line through the balloon position at an angle of o3 Ho 


represents the uppermost altitude of the atmosphere. In addition to 
the notation discussed in I and II, a gas amount matrix U(n, j) is 
introduced here which denotes the gas amount (say in pee acy ee 
along the line of sight n from an 


1 
altitude Hy The monochromatic absorption coefficients, k (i), are 


in layer j, between we and Hi, 


calculated once for each layer using the effective tempterature and 
pressure (as in I and II). 

For the first path from H, at 6. a one-layer (from Hy to Hy) 
monochromatic absorption is assumed according to 


-ky (1)U(1, 1) 
A,(v,H,) =l-e , (6) 


where k (1) is the absorption coefficient at frequency V, accumulated 
from the individual absorbing lines (the first layer). U(1,1) is the 
gas amount in the first layer to be found by the iteration scheme 


(as in I and II), For path 2, U(1,1) is reslanted according to 
U(2,1) = U(1,1) sec 9,/sec 8, ; (7) 


and the absorption in Eq. (6) is recalculated with this new gas amount 
and denoted by A,(v,H)), the absorption in layer 1 along path 2. The 
absorption equation for path 2 is given by 


-ky (2)U(2, 2) ~k,(2)U(2, 2) 
A,(v,H,) = [l-e ] + A,(v,H,)e , (8) 


and the iteration procedure derives U(2,2). Similarly, for path 3 the 
two previous gas amounts are reslanted to derive U(3,3). In general, 
for path n from altitude Ho at | the previous (n-1) gas amounts are 


reslanted according to 
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[eigenen Sn 


| 
| 


= 
en a RR ch a 


U(n, j) = U(j,j) sec @ /sec a for jf = lisscsgmek = (9) 


with which Ay Hep is calculated. The absorption equation for path 
n is 13 


-k (n)U(n, n) ~k (n)U(n, n) 
AL H) = [l-e ] + A, He = (10) 


for n = 2,...,N. 


All Ho » is the absorption of layer n-1 along path n. U(n,n) is 
the gas amount in layer n derived successively. Note that after 
the corrections for the different lines of sight have been made with 


Eq. (9), Eq. (10) involves only two consecutive layers (as in I and II). 
2. Low Sun 


The geometry for low sun data, where both the atmospheric 
curvature and refraction have to be taken into account, is shown (with 
distorted spherical scale) in Figure 7. Here,the balloon is at constant 
float altitude (H,) and each successive low sun path defines a new, 
deeper layer, between two consecutive minimum altitudes. 

The formulation used above for high sun data can be used for 
low sun data provided that 1) the actual airmass distributions along 


the paths are used instead of the sec @ factors, and 2) a single 


effective pressure value can be assigned to each layer (see below). 
Figure 7 shows the first path, which is above the horizon and 

denoted by ray 0, defining the first layer, layer 0, between H, and 

Ho (the top of the atmosphere) with the total airmass along the path 

denoted as M(0,0). This can be used to derive the corresponding 

gas amount U(0,0). The second path, ray 1, passes below the horizon, 

defining the first inner layer from H(-HT™") to Hy (=H observer) 

with the corresponding airmasses 2M(1,1) and M(1,0). The factor of 

2 is included since in the figure only one-half of the path through 
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layer 1 is indicated as M(1,1). U(1,0), the gas amount between H 


0 


and Hy along ray 1, replaces U(0,0) according to 


U(1,0) = U(0,0)M(1,0)/M(0,0) . 


Substituting into 


~k (0)U(1, 0) 
A,(v,H,) = lee , (11) 


and then Ay, H,) into 


-k (1)0(1, 1) -k (1)U(1, 1) 
A,(4H,) = [l-e ] +A,(H,)e 4 (12) 


yields 


U(1,1), the gas amount in the layer H -H, along ray 1. 


1 


Consecutive paths of n@2 each serve to derive a new gas amount 
U(n,n) for the innermost layer n defined by the boundaries So", 


ne (i.e., H to H_), after the outer n gas amounts have been 
n-l ntl n 


adjusted according to 


to 


U(n, j) bs U(j, j)M(n, j)/M(j, j), j=0, eee »n-l, (13) 


This equation applies for n=l1,...,N. Here, U(j,j) and M(j,j) are 
the gas amount and airmass for the previous innermost layers 
1,...,n-1 and M(n,j) is the airmass along path n outside the inner 
layer n. 

The absorption coefficients for a low sun path can be calculated 


as in Eq. (1). The integrations are made along the appropriate sections 
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of the ray path with the ray tracing program. Strictly speaking, when 
a given layer is traversed by different low sun paths (as in Fig. 7), 

a different value of P ote occurs for each path through the layer. How- 
ever, it was found that for the analysis of the type of experiment 
reported here, a single P ate averaged over all of the paths can be used 
(see below). Typical cases are shown in Tables II and III. Without 
this approximation the absorption coefficients for a given layer, which 
depend on P ase’ would have to be recalculated for each path. 

Tables II and III show the airmiss distribution and the effective 
pressure values for two typical cases, at 30 km and at 40 km. The air- 
masses again are given in units of 1 airmass = 2. 153x10°° molecules / 
oni: The pressure is given in Pascals (107 Pascal = 1 mb). The 
upper part of each table is constructed for the case in which the first 
path is below the horizon, thus starting with ray 1 and M(1,1) as 
the total airmass along the first ray. It shows the inner layer bound- 
aries as defined by the minimum altitudes Hain and airmasses and 
effective pressures along the different low sun paths. For rays 
2,3,...N, the airmasses are listed for one-half of the path. A 
column average of the effective pressures is shown in the lower part 
of the table for each of the inner layers. A pressure factor for each 
of the inner layers is defined by the ratio of the average effective 
pressure to the sum of the atmospheric pressures at the boundaries. 

A similar effective pressure and pressure factor are defined for the 
outer layer, H 


to Ho: with the pressure at H, negligible. Also 


shown in the pa part of the tables are siaectints pressure values 
for typical rays just above the horizon. The pressure factor in this 
case is defined as the ratio of effective pressure along the total path 
to the pressure at the balloon altitude. 

It is apparent that the pressure factors depend on the balloon 


altitude and the width of the layers as defined by the low sun paths. 
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A study of this type of calculation showed that, for balloon applications 
and layers that are not thicker than 6 km, a single effective pressure 
value can be determined for each layer which is within +1% to + 5% 
of the individual effective pressures for the different rays. 

These tables are based on a midlatitude summer atmospheric 
profile, but the use of pressure factors permits the application of 
such tables to the pressure profile at the time of the measurement, 
The effective temperatures listed in the tables were calculated in the 
same manner as the effective pressures, but for most practical 
applications, a simple average of the temperatures over the layer 


of interest is sufficient. 
Cc. Photochemically Active Species 


The application of the above methods to photochemically active 
species, such as NO and NO,, is complicated by the time variability 
of the constituent during the sunset period. This time variability 
causes the distribution of the constituent to be asymmetrical with 
respect to H ain’ which results in asymmetrical absorption distributions. 
The effect is obviously largest for the largest solar zenith angles 
and for the most photochemically active species. Thus, this effect 


is quite small for NO,, but very pronounced for NO. In order to 


determine the snaaae Me of the effect on these constituents (which 
are taken to be illustrative of such species), time-dependent profiles 
‘for NO and NO, were used to determine how this variation would 
affect the interpretation of such spectra. The technique used to 
generate these time-dependent profiles is given by L. R. Megill 
of Utah State University in Appendix A. 

For a simplified discussion, we assume again the weak line 
approximation and choose the NO A- doublet [Goldman and Schmidt, 


1975] at 1900.08 cm7! with S = 1.0x107!? em7!/(molecule cm”) at 


300°K and ground state energy of E"' = 80.3 em", and the NO, four 
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lines group [Goldman et al., 1975], including spin-splitting, at 1604, 57 


19 cm Meisleods wee at 296°K with E" = 


em! with S = 4, 6x10" 

97.8 aunt: These are typical spectral features used for high altitude 

observations of these trace gases. The temperature dependence of the 

line intensities is included according to standard equations. Figures 

8 and 9 show calculated integrated absorption distributions for NO 

and NO, at 93° from 30 km. Within the (small) temperature dependence 

these distributions are proportional to the gas amount distributions. 
Figure 8 shows that most of the NO absorption is distributed 

widely over the portion of the path between A nin and 40 km in the 

sun direction. If the data are analyzed on the basis of a non-time- 

dependent distribution using the full analysis presented above, the 

absorption in a given layer would be attributed to the portions of the 

path on the sun side and on the platform side. In some cases there 

is no absorption on the "dark" side, which ca. result in a 50% error 

in mixing ratio derived for the layer. If one of the approximations is 

used (e.g., a one-layer calculation which assumes that the absorption 


takes place in a narrow layer above H nin) the inferred NO mixing 


ratio could be in error by an order of magnitude. In the case of NO,, 


as seen in Figure 9, the distribution along the path is not so asymmetrical, 


so that the error in the case of a full analysis is much less. However, 


as mentioned above, approximate calculations can result in large errors. 


IV. Personnel 


The major calculations performed in this report represent exten- 
sions of programs developed by the Upper Atmospheric Physics Group 
at the University of Denver. The major effort in the development of 
these programs has been by A.Goldman. He has also been responsible 
for the extensions which address the problem of the interpretation of 
sunset spectra. He has been assisted in this effort by G.R. Cook, D.K, 
Rolens and R.S.Saunders. The figures have been prepared by C, M. Bauer. 
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V. Reports and Publications 


As indicated in the text, this report contains summaries of three 


reports issued as part of this project. These are: 


Fe 


Goldman, A., G.R. Cook, D.G. Rolens, and D.G. Murcray, On 
the Interpretation of Infrared Sunset Spectra for Altitude Distri- 
bution of Atmospheric Trace Constituents, Scientific Report 

No. 1, Contract No. DOT-FA77WA-3949, Federal Aviation 
Administration, Dept. of Transportation, by Department of 
Physics, University of Denver, Denver, Colorado 80208, Nov. 
1977. 


Goldman, A., G.R.Cook, D.G. Rolens, D.G. Murcray, L. R. Megill 
and W. Abdou, On the Interpretation of Infrared Sunset Spectra 
for Altitude Distribution of Atmospheric Trace Constituents, by 
Department of Physics, University of Denver and Center for 
Research Aeronomy, Utah State University, 1977. 


Goldman,A., and R.S.Saunders, Analysis of Atmospheric Infrared 
Spectra for Altitude Distribution of Atmospheric Trace Consti- 
tuents, I. Method of Analysis, in press, J.Q.S.R.T., 1978. 
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Appendix A 
"Temporal Variations of Minor Constituents at Twilight 
L. R. Megill 
A one-dimensional, time dependent model is used in the present work 
to describe variation in the distribution of a number of important 
gaseous compounds (e.g., NO, , C10, , HO, , and CFM.) from 10 km above the 
ground level up to 80 km, during sunset. In the calculations, the solar 
zenith angle varies from 80° to 95°. In the model used the continuity 
equations describing the rate of change with time of the concentration 


of a species x is given by the equation 


ee (1) 


where P and L are the photochemical production and loss terms, respec- | 
tively. The vertical transport processes are ignored in the present 
work since the 15° solar depression considered by the calculations 
corresponds only to about 2 to 3 hours during which integration of 
equation 1 is performed. This time is short compared to the time con- | 
stant for transport effects to be significant. 


The accuracy of the results depends, among other things, on the 


aren RAO 


number of chemical and photochemical reactions used in the calculations. 
The data present in literature, however, limit the number of reactions 
used. Figure 1 represents the basic reaction schemes considered in the 
present study and a listing of the reactions involved is renresented in 
Table 1. The rate constants of the chemical reactions used are mostly 
recommended in NASA report number 

' Other data collected for the completion of this work include: 
1. A set of starting profiles for each species involved. 


2. The photodissociation rates for the photoreactions. 
A-l 


The starting profiles for most of the species were found in the lit- 
erature only up to 40 to 60 km. These profiles were extrapolated up to 80 


km using the scale height equation ; 


e WH 


N=N, (2) | 


where H is the scale height, N, is the concentration at altitude z and No» 


the concentration at a height h above the altitude z. 


Those starting profiles, especially of the sunset sensitive 
species, are only approximate since they do not represent the actual F 
profiles at the starting time of the calculations. 
Sets of starting profiles which represent greater internal consistency 
could be obtained from the available ones by successively integrating 
the continuity equations for 24 hour periods until the essential diurnal 
effects are properly reproduced. This however, has not been done in the 
present work and is intended to be accomplished in the next phase of this 
work. 
The photo-dissociation rates were calculated by a computer code written 


by Harris and Adams. The ca:culations are based on the equation where 


ho - 
5,0,A,¥,2) =f 9, Q) 6 Me “nm Matin (3) 
Ay 


where i denotes the species, m denotes the major absorbing species, jj 
is the dissociation rate per molecule of species i (in sec”), xy and 


A, are the wavelengths bounding the region of interest, x is the zenith 


angle of the sun, z is the altitude, 9; (o,) is the photon absorption 


cross section of species i(m), ¢ is the solar flux at the top of the 
atmosphere and NOX 2) is the column density of species m along the path 
from altitude z to the sun at zenith angle x, i.e., Na = LIN, Secy dz. 
Single scattering is the only effect considered by the code in calcula- 
ting the photodissociation rates. 0, and 0; are considered as the major 
absorbing species. Their density profiles, were taken respectively from 
U.S. Standard Atmosphere (1962) and from Adams and Megill (1970). The 
absorption cross sections of 0. and Oz are represented graphically in 
Figure 2 and 3, respectively. For consistency in the calculations the 
density profiles and the absorption cross sections of 0 and Oz used in 


the present work are the same as those used by the code calculating the 


photodissociation rates. 


Calculations 

The coupled differential equations represented by Equation 1, were 
solved by a computer code written by Adams and Megill (1970). A modified 
version of this program is used in the present work. The integration 
routine of the code is based on a technique which uses an exponential vari- 


ation during the time step. The accuracy of a step is checked by backward 


integration. Agreement to 1 part in 10° is required by the program. The 


step size is automatically adjusted as needed to keep this accuracy. An 
initial step size is taken as 0.1 sec. This could be halved or doubled 
several times to: maintain the required accuracy. 

The output of the calculations gives the distribution of the species 


at a number of points representing space and time. These points constitute 
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a 6 x 151 gridpoints representing height variation of 2 km (ranging from 


10 km to 40 kn) and solar zenith angle variation of 0.1 ranging from 80° 


to 95°. The geometry of the grid is sketched in Figure 6. This grid 


allows the representation of isodensity plots for any of the species in- 
volved. Those for NO, NO, C10 and Cl are represented in Figure 4,5,6, 
and 7, respectively. 


TABLE (1) 


a - Chemical Reactions Set 


| | __ ae _—-—  e Costa | Reference 
1. NO + 0, -----—— NO, + 0, 2.1 x 10° exp - 420 | | 
2. NO, + 0 ——-—--+ NO, + 0, 9.1 x 107? | 
3. N+0, —-——+NO+0 5.5 x 101? exp -3220/T | 
4. N+NO —-—-—+N, +0 8.2 x 101) exp -410/7 | 
5. NO, + OH +M —-—-—HINO, + M~—|See below (*) 
6. OH + HND, -—-—---» HO +NO, =f x 1074 | 
1. N + NO, ——--—-—+ N,0 +0 2.0 x 107+ exp 20? | 
| 8. NO+O0+M —10, +M 1.55 x 107 ep S84 | 
| 9. NO, + HO, +.Ny ———-—+ HOQNOs#N2 [1.0 x 1072? | 
10. N + 0, --—--—---—> NO + 0, 5 x 107? exp “890 | 
| 1. NO, +0; ——-—-—+ NO, +0, [2.2 x 10719 exp 2450 ‘ 
12. olp + N,0 —----—— N, + 0, 5.5 x 104 : 
13. olp + N,0 -—-—--— NO + NO 5.5 x 10 tt g 
| 14. olp + H,0 ----——> OH + OH 2.3 x 102° 
15. ol + CH, —--——> OH + CH 1.3 x 101° ; 
16. Ol) +No +M-—-——— NO +M = [3.5 x 10727 | 
17. oly + Hy ——-—-——> OH +H 9.9 x 10+ | 
| 18. oly +N, ———--— 0 +N, 2x10! ep HY | 
19. olp + 0, -——--—-—— 0 + 05 2.9 x 107+ exp 67/T | 
20. olp + 0, ——-—-—> 0, + 0, 1.2 x 102° 
| 21. oly + 0, ——--——-+ 0, +0 +0 {1.2.x 107 
) 22. oD + Hel -—-----— OH + c2 1.4 x 10720 | 
! 23. olp + CFel3 ——-——-—+ Cel telo [2.3 x 107° 
24. Olp + ely ——-—-—+ CF el +10 [2.0 x 107° 
25. OH + HO, —————+ Ho +0, = [3x 107 
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TABLE (1) continued 
a - Chemical Reactions Set 
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Nie) Reaction Rate Constant Reference 
26. HO, + HO, ——-—--—+ H,0, +0, | 5 x 107 exp 200 
27. OH + OH -----—--+ 1,0 + 0 1.0 x 107 exp 9/7 | 
28. HO, + 0, -——--——+ OH +. 04.0, | 1.0 x 1077? exp 3% | 
29. OH +0, ——--——> HO, +0, 1.5 x 107)? exp =2000 i 
30. OH +0 —-——---—+ H + 0, 1.0 x 10729 exp ~250/p 
31. O + HO, -——---——+ OH + 0, 1.0 x 10729 exp 7259/7 | 
32. 0 + H,0, ——-—-—> OH + HO, 2.75 x 10)? exp ~2125 7 | 
33. H +0, +M—-—-+ HO, +M 2.1 x 1073? exp £220 
34. NO + HO, -----——> NO, + OH 1.5 x 107'* exp HOOT 
35. H + 0, —-———-+ OH + 0, 1.2 x 107"? exp 82 | 
36. OH + OH + M —--——> H,0, +M 1.25 x 10-7? exp #200 ‘ 
37. OH + HO, —— + HO + HO, 1.0 x 10 2+ exp 7/977 B 
38. OH + CO —---——+ 00, +H 1.4x 1073 E 
39. OH + CH, ——----—> HO +.cH, | 2.35 x 107)? exp “10 3 
40. el +03 -——--—-—+ clo +0, 2.7.x 101) exp 27/7 & 
41. elo +0 --——-—> el + 0, 7.7% 10 exp 7130/7 
42. elo + NO -—--——-> el + NO, 2.2x 10 1 | 
43. OH + Hol----—— ——+ cl + 1,0 3x 107)? exp ~425 7 | 
44. el + CH, ------- — Hel + CH, 7.3 x 101 exp ~1260/- | 
45. el +O, ——--—+ el +0, $230 
46. elo + NOgtN,—--—+ cloNO, +Ny | 5.1 x 10°? exp 133 | 
47. © + Big) eeeeeees 9014 OH 1.4.x 10722 exp 3370/7 | 
48. ol + Oh ——--——» 0+ Het 1.0 x 1071} exp =2370 | 
49. cl +H, ———————+ Hel + 3.5.x 1072) exp ~2290 7 | 
50. el +H,0, -——-—+ Hel +0, | 1.7.x 107)? ep 4/1 | 


TABLE (1) continued 


a — Chemical Reactions Set 


Rate Constant 


el + HNO, -—-—--— Hel + N03 1x aa exp “2870 pp 


A 
03 +0 <---> 05 + 05 1.9 x 197} exp -2300/,, A 
0, +O+M——-—-—> 0, +M 1.07 x 107 *4exp*320 


0+0+M ———— 0, +M 8.0 x 107° 


Hy oo ——__—_ > f+ H,O 2.3°% 19721, iia 


-11 
HO 2. Hy + 05 1.0 x 10 


=31 


Olp + BH -21g + 8 éixie 


-31 
HO>+NO5+M ~~~ ———+H05 +NO5 1M 1.97 x 10 


eS HO + NO 3x10 exp 


-390 
t 


2 
=9 = 
HO,NO,+e2 -—-----—> Hel + 0,+NO, 1.7 x 10 “exp 


oO 


HO,NO, + 0 ---—-—> OH40,4NO, 2.75x10 1? exp 
CO + CH -~——---—> H + ©, 
CH,0 + 0, ---—--—> CH,0 + HO, 

5.3x 10. 


N+ HO, -~-—--—-—+ NO + CH 2.0 x PT Bors 


HO,NO,, eae ecisrmint cacti H,0 + 0, 1.0 x 10 


NO, + NO, +M—--—> NjO, + M See below (**) 
CH,O, + NO --—--—+ CHO + NO, 1.5 x 10 exp 9p 
O + OH ------— 1,0 + al 1.4x 107} 


2 
CHO +O - 


11 


on OP FF Pp FP YP GD 


CH 


> YP PY 


os HO, + Co 1.7 x10. 
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(A) Crutzen, 1977. 
(B) Atmospheric Environment, Volume 9, 1975, page 1045. 


* k for reaction 5 is given as 


k 
[e) 
pa ae 


k, and k, are the low pressure and the high pressure limits 
for the rates, respectively. 

[u] is the air density 

k, = 6.5 x 10°77 Fanos} 


k, = 1.0 x 107 ans) 


**K for reaction 69 is given by the expression 


pe ky 3+ vu] uae 
I+ svn] 3+vu i a 
k, 
Peng Ue 
k, = 1.9 x 10%cm°s™* 


k, = 3.8 x 10 eens t 
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b - The Photoreaction Set 


a ee ee 


02 


0, + hv ----—-—-+ 


+ hv ee ee mee 


Coe ae 
NO +hvy -—------—> 
Pt ee eee 
N50, + [hh 
py thy 
H NO, + hv —— 
i he eer 
elo NOj+ hv --—--—> 
oe ft ih 
Cel, + hy <r 
Gets + hy 
CFycl, + hy --—--—> 
CieZ5F + hy <---> 
CHc2 + hy ---—-—-> 


CH,02H + hy ----—-> 


HO,NO, + hy ---———> 


N,0 + hy —-—--——> 
oot hy = 


Hel + hy --—-———+ 


0+0 
O+ 0, 
Oy, + 02 
NO +0 
NO, +0 


NO, + NO. 


2 3 


OH + OH 
OH + NO, 
OH +O 
NO, + clo 
el +el 
eelg +cl + el 
el + CFel., 
et + CFocl 
CHelF + cl 
CH, + cl 
CH, O + OH 
N, + Oy 
el +o 


H+ el 


J values for H50. are used 


J values for H,0, are used 
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APPENDIX B 


On the Interpretation of Infrared Sunset Spectra 


for Altitude Distribution of Atmospheric Trace Constituents 
by 


A. Goldman, G.R. Cook, D.K. Rolens” and D.G. Murcray 
Department of Physics 
University of Denver 
Denver, Colorado 80208 


Scientific Report #1 
Contract No. DOT-FA77WA-3949 
Federal Aviation Administration 
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Abstract 


Quantitative interpretation of sunset spectra in terms of mixing 
ratio altitude profiles of atmospheric trace constituents requires, 
amont other things, a detailed consideration of the air mass distribu- 
tion along the sunset ray. Air mass distribution tables are generated 
for such infrared sunset spectra by a ray-tracing technique. The ap- 
plication of these distributions to the interpretation of sunset spectra 
is discussed for both sunset insensitive gases and sunset sensitive 


gases such as NO and NO,. 


November 1977 


* 
Laboratory of Applied Mechanics, Denver Research Institute, Denver, 
Colorado 80208 
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Le Introduction 

A powerful infrared technique for the detection of trace constituents 
is to observe atmospheric transmission spectra from an airborne plat- 
form using the setting sun as a source. The main advantage of this 
technique is the long path achieved through the atmosphere below the 
observation point, which increases the spectral sensitivity of the method 
to atmospheric trace constituents, In fact, within the last decade, 
several trace constituents have been measured in the lower stratosphere 


NO,, NO, CF,Cl,, CFCl,, CCl,, HCl, 


by this technique, e.g. HNO gute: 


and HF. 


3’ 3’ 4’ 


Quantitative derivation of mixing ratio altitude profiles of trace 
constituents from sunset spectra requires a detailed consideration not 
only of the spectral absorption lines but also of the geometry of the 
optical path through the atmosphere which, in turn, is largely affected 
by atmospheric refraction. In a low sun experiment the path goes through 
a minimum altitude, around which a significant portion of the total air 
mass is concentrated. Thus, to a first approximation it has been 
customary to assign the observed absorption to the vicinity of the mini- 
mum altitude. In the present paper, we study in detail the air mass 
distribution along such a path for an isotropic spherical atmosphere 
divided into thin homogeneous shells. These distributions are then 
applied to the interpretation of sunset spectra for sunset insensitive 
‘gases and also for sunset sensitive gases such as NO and NO,. 

The air mass distributions for the layers of an isotropic atmosphere 
are symmetrical with respect to the minimum altitude and have 
significant widths which depend on the zenith angle and the altitude. 

In contrast, current time-dependent photochemical models predict that 
the NO and NO, distributions are asymmetrical with respect to the 
minimum altitude, especially NO as its concentration can change drastically 
during the sunset period. 
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The tables of air mass distributions produced in the present study allow 

an accurate quantitative layer by layer derivation of the mixing ratio profiles, 
assuming sunset insensitive gases. Common approximations are 

discussed and compared with the exact method for the case of the 

weak non-overlapping line approximation, The same method can be 

extended to a line by line - layer by layer calculation, The air mass 
distribution tables are then used to examine the difficulties associated 

with the interpretation of sunset spectra of sunset sensitive gases, 

It will be shown that a derivation of an NO profile can be quite 


meaningless, while a meaningful profile for NO, can be derived. 


2 
2. Air mass distributions 

The air mass distribution along a sunset path is derived from 
a ray-tracing computer program, including refraction, described 
in detail in Snider and Goldman (1975) and Snider (1975). A standard 
midlatitude summer atmosphere, as tabulated by McClatchey et al. 
(1972), and a wavelength of 5 ym have been assumed for the present 
work. Typical air mass distributions, as well as the accumulated 
distributions, for an observer at 30 km and apparent zenith angles of 
92° and 94° are shown in Figs. 1-4. Table 1 is an example of the 
corresponding tables and presents the numerical values for 94°, 
Similarly, air mass distribution tables and figures were made at 
0. 5° steps from 90. 5° on, until the ray strikes the earth's disk, for 
observer altitudes 15, 20, 25, 30, 35, 40, 45 and 50 km. These 
are available along with NO and NO, absorptions in the Appendix C. The 
0.5° step is dictated by the sun's disk (32 min of arc) and typical 
spectral scan times (~1 minute). The computer programs used 
can easily be rerun for a different set of altitudes and sun angles. 

For a given observer altitude (in kilometers) and an apparent 
zenith angle (in degress), the air mass distributions presented in the 
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tables and figures are divided into two parts. The first part 

includes the values from the minimum altitude, A nin’ to the observer 
altitude; the second extends from ae towards the sun,up to 100 km. 
These two parts are separated by a vertical line at A nin in the figures 
and by a space in the tables. The total air mass listed on each 

figure is along the total path which includes the standard atmosphere 

up to 100 km altitude. One air mass is 2.153 x 107° molecule /cem*” 
(Snider and Goldman, 1975). The ray-tracing computations have been 
performed in 250 meter steps for the 0 - 40 km altitude range, 

1 km steps for 40-70 km and 5 km steps for 70 - 100 km, as explained 
by Snider and Goldman (1975). The figures and the tables have been 
calculated in 1 km atmospheric shells, so starting at Fain the first 
step can be less than or equal to 1 km until the nearest 1 km shell 
boundary is reached. Thereafter, the steps are the 1 km shell 

steps. Figure 5 is an example of ray geometry and air mass 
distributions along various paths. It is seen in Figs. 1-5 that the 

air mass distributions are symmetrical around Hain and have signifi- 
cant width. There is always less than 50% of the air mass between 
Hain and the observer and more than 50% between Hain and the sun. 
The air mass tables are extended to an arbitrary altitude, H 


f’ 
is 10 km above the altitude of the observer in the sun direction, so 


which 


that the residual air mass is quite small (but not always negligible). 
The truncated air mass listed in Figs. 1-4 is the air mass from Hain 
to Hy and the residual air mass above H, can be determined by 

a comparison with the total air mass listed. The amplitude in 

the air mass distribution plots (Figs. 1 and 3) are relative to the 

total air mass. Figures 2 and 4 show accumulated air mass distributions 
as a running accumulation of the 1 km step values on both sides of 


Henin which are also relative to the total air mass. It should be 
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noted that each tabulated value in the tables represents the 1 km shell 


with the printed altitude as its lower boundary, Thus, the last line 
tabulated in Table 1 as 39.0 km represents the 39 - 40 km shell, 

The tabulated zenith angles are along the path, at the intersection 
points of the altitude boundary listed. An example of these values 
along the 92° path are shown on Fig. 5. Thus, A nin is always at 


90° and the angles in either direction from Fini complement one 


in 
another to 180°, one increasing and one decreasing from the 90° 
value. The tabulated pressure (in Pascals, 10° Pascal = 1 mbar) and 
temperature are average values over the shells. The effective 


temperature and pressure in the tables are defined by 
P ott = J Pdm/ y dm and Tore = J Tdm / J dm , 


where dm is the air mass increment, P the pressure and T the tem- 
perature, as explained by Snider and Goldman (1975). The integrands 
for these expressions, after the small step computations have been 
accumulated over 1 km steps, are given in the tables under the columns 
A(P-dm) and A(T-dm). The air mass increments are listed under AM. 
The 1 km step values are accumulated along the path from Aynin in both 
directions under the column ACCUM &M. These are the values shown 
in Figs. 2 and 4 as the accumulated air mass distributions. The geo- 
metrical path increment is listed under A(PATH) and is accumulated 
along the path under the ACCUM A(PATH) column. Thus, for 30 km 
at 93°, the geometrical path from the observer to Hain is 338 km and 
for 94° is 457 km. It should be noted that with such distances the 


isotropic atmosphere approximation can be quite poor if the path 


reaches the troposphere above heavily polluted cities. 


PT LL 
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For a uniformly mixed gas in an isotropic atmosphere, the gas 
amount distribution is proportional to the air mass distribution. 
Further, in the weak-line approximation the gas amount distribution 
is also proportional to the integrated absorption distribution, within 
the temperature dependence of the line intensities. Thus, the inte- 
grated absorption distributions of such gases have approximately 
the same shape as the air mass distributions of the type displayed 
in Figs. 1-4. Significant deviations from these shapes indicate 
altitude dependent mixing ratios if the shapes remain symmetrical 


and nonisotropical gas distributions if the shapes are asymmetrical. 


3. The interpretation of infrared sunset spectra 
Denote the observed integrated absorption at different apparent 
zenith angles z ese, BY A A 


z » for a given observer 


ae ie Sagi 


altitude, h Without a loss of generality, assume that hy is an 


integer Oe altitude. For simplicity of the discussion, the weak 
line approximation is assumed, so that the integrated absorption 
A, is given by ‘Atces = § tom" /molecdls on} Ginclecaiej cus), 
where S is the integrated line intensity and U is the absorber amount. 
The air mass distribution tables give atmospheric amounts U as MU 
where M is the tabulated air mass, and Ue = 1 air mass = 2,153 x 107° 
molecule fens”. Trace constituent amounts will be described by BU 
where 8 is the gas volume mixing ratio. 

The tables can be used for a layer by layer derivation of a 
minor constituent altitude profile in the following way. The atmosphere 
is divided into even one km shells, although Fig. 5 shows 2 km shells 
for simplicity. The tables show 1 km shells which can be accumulated 
over two or more km as required. Starting with the first path 
and z 7 90° (required to observe the weak absorption of the trace 
constituent), a constant mixing ratio, By» can be derived for the 


1 
* ) of the path according to 


B-6 


atmosphere above H_., 
mi 
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where M, is the total air mass along path 1, This B, will be considered 


1 
as the mixing ratio for layer 1, between h, and h,, as well as for the 


residual atmosphere above h, in the sun direction, It should be 
noted that A nain can obviously fall between the even km boundaries, so 


that h, will be the same altitude as A nin or the next nearest even 


km below it. Subsequent inner layer boundaries are defined by 
the minimum altitudes of the corresponding zenith angles, Z.» S80 that a 
layer can include more than 1 shell with increasing zenith angle. For 
example, with the observer at 30 km, the inner layer defined by z = 92° 
and z = 93° includes five 1 km shells, This is shown as Layer 3 
on Fig. 5. 

A mixing ratio altitude profile is then derived by the successive 


equations: 


n-1 


Pg = ee a RMbiar 74 
2 ag eee : n=2, 3... thi 
sU_M 
o on,ntl 


Here, Bo is the constant mixing ratio for the layer defined by H 


(n-1) (n) , 
to hoy the even km boundary below A nin . oe et] 


portion of the air mass between altitudes hy and bea for path n. The 


upper index n denotes the different paths for a given observer altitude. 


min 


The M,. ee values are those tabulated in the air mass distribution 


tables every 1 km accumulated over the proper layers. In Eq. 2 the 


absorptions in the outer layers, each one of which is crossed twice. 


A smaller correction is a = psu mM '") » where M. is the residual 
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summation in the numerator represents the correction to A, due to the | 
J 


RE crane 


| 
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: 
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Seal ina altel aaa octal atin casteN ates air ons t oo Snes Nene aiaigilis, alt ci cit ly nits a 


air mass outside h, in the sun direction, The inner layer n is 


1 
crossed once, but note that its air mass is twice the air mass from 


: to h 
min n- 


\ 
The corrections to AL are by no means negligible. However, a 


commonly used approximation is to neglect these corrections in the 
(n) 

n, n+l 

total airmass along path n, and assign the resulting mixing ratio, 


numerator and replace M in the denominator by M_; the 


6, ~ 30M (3) 
on 


to H ., (n) 
min 


motivated by the fact that a significant portion of the total air mass 


or a small altitude range above it. This approximation is 


along the path is indeed concentrated near Ain’ It can be a good 


approximation only if the effect of neglected corrections in the numerator 


; (n) 
are compensated for by the ratio M, / ~~ as so that 
> (2) ,(a))__n 
(A, - se 2 Ae Mi, tot - 4 | ener a 
- n, n+l 


Examination of the air mass distribution tables shows that for many 
practical cases, this approximation can lead to errors from 5% to a 
factor of 5, as will be demonstrated below. 

Another commonly used approximation is to consider the absorption 
, and the corresponding total air mass difference, 


-1 
, so that BR. is approximated by. 


difference, A_ - A 
n n 


M_ -M 
n- 


n 1 


ms ; “ee 


Pn OM, = 50, ss 


(n) the 
and assigned to H_. 
min 


as inaprevious approximation. This approximation 


can be quite good if both AL in and M, >M but in many 


2) | n-1’ 


practical cases A, =A and M. is not very large compared to 


n-l 
M 1? unless z2 94°, A combination of these conditions can result 


Py very large error. 

Both of these approximations are a one-layer calculation. Strictly 
speaking, the resulting mixing ratio should then be assigned not to 
A nin ,» but rather to a wider layer, the width of which increases 
as z increases, as is shown by the air mass distributions. If an 
arbitrary definition of a 90% distribution width is adopted, the 
width can be quite large, as is seen in Figs. 2 and 4. Thus, for 
30 km and 90% of total air mass, the width is ~5 km for 92° and 
~7 km for 94°, The exact Eq. (2) allows a finer altitude resolution, 
with the overlapping betw ir mass distributions at di 
zenith angles properly accounted for. A typical altitude resolution 
obtainable by using the entire sun’s disk from 30 km altitude is 
achieved by starting at 90. 5° and proceeding in 0. 5° steps. The 
minimum altitudes for 90.5, 91.0, 91.5, 92.0, 92.5, 93.0, 93.5, 
94.0 and 94.5 are (Snider and Goldman, 1975) 29.7, 29.0, 27.8, 26.1, 
23.9, 21.1, 17.9, 14.1 and 9.7 km, thus defining 9 layers with 
widths of 0.7, 1.2, 1.7, 2.2, 2.8, 3.2, 3.8 and 4.4 km. 

As an example, for the application of Eqs. (1) and (2) and 


the approximations (3) and (4), consider an observer at h, = 30 km, 


1 
and three integrated absorption measurements at z= 92°, z, = 93° and 
z= 94°, This can be visualized with the help of Fig. 5, but the 91° 
path does not pertain to this example and the 94° path is not shown 

(for simplicity). From our tables we find the minimum altitudes of 
H_, oe) = 26.08 km, H_. (2) = 21.15 km andH , (3) = 14.09 km; 

min min min 
and total air masses of M, = 1,581, M, = 3.756 and M, = 12.088. 
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Assuming measured integrated absorptions of A, = ge A, = 
and A, = 2 a which indicate a trace gas ere in a 
layer in the atmosphere, and assuming S = 1 x 10° cm/molecule, we 
find by Eq. (1) that By = 2.94 x 10°78 vol/vol. 

The approximations (3) and (4) use the same B, value. The exact 


B, for the next layer (21 to 26 km) is given by 


A, - am, ,'e,su, - Mg, su, : 
: x = 6.89 x 10°” vol/vol , 


B = 
2 (2) 
SUM), 


whe re mf = 3.08 (between 21,15 and 26 km), 2M,, = 0.480 
(between 30 and 26 km) and m '?) = 0.196. The increased mixing 
ratio in this second layer shows a gas layered in the atmosphere, 
since A, is relatively small. 

Approximation (3) yields B, z A,/SU,M, = 6.18 x i0°8 vol/vol, 
which is 10% lower than the exact B,. This approximation can 
improve significantly for smaller A, values, i.e., for a gas which 


is less layered in the atmosphere. Approximation (4) yields 


A,-A, 


-8 
B, = (M, M,)SU, = 8.54x10 vol/vol, 


which is 24% higher than the exact B,- This approximation will 


be improved for larger A, values. The exact B, for the next layer 


(14 to 21 km) is given by 


(3) (3) (3) 
rhe A, 4 2B) Mis au. bs 2B, M2, ou, : M, A au, h. 
3 su,M,,'” 


5.68 x 197 1° vol/vol , 


(3) (between 26 and 30 km) = 0, 334, am, ,'*) (between 21 
and 26km) = 1.050, i (between 14,09 and 2lkm) = 10. 544 and 
a = 0.160. Approximation (3) gives e. = A,/SU\M, = 7.68 x ie, 


more than 10 times higher than the exact B.. Again, this approximation 


whe re 2M, 2 


is improved for smaller A, values. Approximation (4) for B, is not 


applicable here since A, - A, <0. It is thus seen that the common 

approximations can lead to a very large range of errors as stated 

earlier. The numerical values of this example are summarized in Table 1. 
The layer by layer method described above can be extended 

to a more general line by line - layer by layer calculation over 


a spectral interval of non-weak overlapping lines. This can be done 


similarly to the presentation by Goldman et al. (1973). For the 
first path 


AMy 1 = wn, (5) 


and for the successive inner paths 


n-l n 


AM yn = Mym+2 FS eMym no Myre e™ , 
m=1 k=m+1 . 
WW SZ yisg tesa OG (6) 


where the upper index n denotes the path n, with "infinite" resolution 
spectral transmitivity, oy, k), and absortivity ) (V,k) at wavenumber 
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V of a layer defined by the altitudes hy to Bet" AM, n) is the 
spectral absorption along the total path n. € ‘(UV,n) is the emissivity of 
the inner layer whose mixing ratio B. is being determined, am is 
the emissivity of the residual of the atmosphere outside h, in the 

sun direction, By fitting the observed integrated absorptions to the 
theoretical ones as done by Goldman et al. (1973), a mixing ratio 
altitude profile can be derived. In the weak line approximation 


for a single line, Eqs, (5) and (6) yield Eqs. (1) and (2). 


4, Sunset sensitive gases 


The application of the above methods for sunset sensitive gases 


such as NO and NO, is complicated by the uneven gas distribution during 


the sunset period, These distributions are asymmetrical with respect 
to FH nin’ resulting in asymmetrical absorption distributions. The 
degree of asymmetry increases for larger zenith angles, It is 

quite small for NO,» but very pronounced for NO as will be described 
below. 

For a simplified discussion, we assume again the weak line 
approximation and choose the NO A- doublet (Goldman and Schmidt, 
1975) at 1900.08 con” with 5 102 10°°" om” */(molecdle em} 
at 300 K and ground state energy of E"' = 80.3 em” *, and the 
NO, four Sesiiae ait: (Goldman et a. weg including oe 
at 1604.57 cm § with S = 4.6 x 10 cm /(molecule cm ‘) at 296 K 
with E" = 97,8 cm™/, ‘These are typical spectral features used 
for high altitude observations of these trace gases. The temperature 
dependence of the line intensities is included according to standard 
equations. Figs. 6 and 7 show calculated integrated absorption 
distributions for NO and NO, at 93° from 30 km, Within the (small) 


temperature dependence, these distributions are proportional to the 


gas amount distributions, The latter were taken from a time-dependent 
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photochemical model tables of NO and NO, number densities at various alti- 
tudes and zenith angles, generated by R. Megill and W. Abdou. Additional 
NO and NO, figures and tables for various cases are in the AppendixC. 
Figure 6 shows that most of the NO absorption is distributed widely 
over the portion of the path between Hain and 40 km in the sun direction, 
Thus, if the observed absorptions during sunset are used to derive mixing 
ratio distribution by means of Eqs. (1) and (2), the NO absorptions are then 
assumed to be symmetrical around Hain’ which will lead to a large error 
because there is more than an order of magnitude diffe rence in the NO mixing 
ratios in either direction. Here it should be noted that a given shell has dif- 
ferent mixing ratios for different paths due to the time dependence. As an 
example of such an inversion of the data, we consider again the case of 30 km 
at 92", 93° and 94°, The pre-calculated mixing ratios from the photochemical 
intice’ for the three layers of 30-26 km, 26-21 km and 21-14 km (as de- 
1, Bat )=7.95x10" 
=1,39x 


fined by the given zenith angles), are gM 6.09x 10° 
and g (12 22x10 ae for the H_., to observer direction, and B ay 
3 -i0 min a? 1 
10°”, g,.=4.54x10° and g,{!)= 1.90x10°°" for the H_ to sun direction. 
2 3 . min -9 
The corresponding inverted g values, by Eqs, (1) and (2), are B,7}: 13x10 /, 


R= 2.55x 10° a and B37: 68x10 sah It is apparent that while the averages 
of 2 'D) (II) 
B; 


and B. are within 5% and 50% from the inverted Be the variation 
between B; and B; 


(I)" (II) 


can exceed an order of magnitude (depending on the 
photochemical model used). Such possible variations within the shell are not 


accounted for in our inversion method and the resulting mixing ratios are not 


always meaningful. Further errors can be introduced if Eq. (3) or (4) are used. 


In particular, 


* For simplicity the integrated absorptions were calculated without the tem- 
perature dependence (T) of the line intensity S(T) so that A, = 3, 842x 1073, 
A,=3. 320x Ts and A ,=3. 280 x 10" "ein" % These A values differ somewhat 
from the tabulated values in Appendix c which do include the temperature 
dependence of S. The temperature independent values have been chosen to 


simplify the example of using Eq. (1)-(2) where, strictly, S should vary with 


each layer according to temperature, 
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the precalculated integrated absorptions for NO along the 92°, 93° 
94° paths are 0.00384, 0.00332 and 0.00328 in” which are almost 
equal if an experimental error is allowed. 

The NO, distributions are only slightly asymmetrical around H nin 
as demonstrated in Fig. 7, so that the inversion process will be much 
more accurate. In fact, for the same example as above of 30 km and 
92°, 93° and 94°, the precalculated integrated absorption values of 


0.0765, 0.101 and 0.124 eae (again calculated without the temperature 


aps aig of the line intensities), we find B, = 5.82x 1077, ae = 
2.85x 10°” and 3 (1). 1,04%10°? for the H i, t0 observer side, and 
B') 25.05x10%, Bo) = 248x107 and =, = 8.51 x 1071 
Be the H_. to sun day The inverted values are .* 4.89x 10° me, 


B= 2.25% 10° and A= 7.35% 107°: Sith Witéws "a mieacingtat 


inversion for NO, despite the slight asymmetry. See Table 1 for a 


2 
summary of these comparisons. It may seem odd here that the B values 
I II 
are smaller than both p ) and p! ) but this is due to the fact that 


B (I) 


and gill) were obtained by averaging the mixing ratios of 1 km 
shells over each layer. The inversion is limited in altitude resolution 
and the layers include several 1 km shells. 

It is thus seen, that a time-dependent gas distribution model can 
significantly affect the quantitative interpretation of sunset spectra. 
Since the sunset long path is needed in many cases to detect the other- 
wise very weak absorption of a trace constituent, a time-dependent 
model must be a part of the data interpretation. The layer by layer 
analysis can be used to test such models, even though the results will 
not always be unique. High sun or constant low sun observations should 


then be made to eliminate time-dependence effects. 


CF Summary 
Air mass distributions along sunset paths have been generated 
and used to exarnine some problems concerning the interpretation of 
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sunset spectra in terms of mixing ratio altitude profiles of atmospheric 
trace constituents. It has been shown that the distributions have to be 
taken into account for both sunset insensitive and sunset sensitive gases. 
In the first case, commonly used approximations which neglect 
these distributions can lead to significant errors in the derivation of 
mixing ratio altitude profiles. The second case is further complicated 
by the asymmetry of the gas distribution with respect to the minimum 
altitude of the path. The tables in Appendix C provide a base for 
correcting the interpretation of the infrared spectra of such sunset 


sensitive gases. 
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Table 1. 


Table 2. 


Figure and Table Legends 


Air mass distribution for | km shells along a 92° path from 
30 km. 


Accumulated air mass distribution over 1 km shells along a 
92° path from 30 km. 


Air mass distribution for 1 km shells along a 94° path from 
30 km. 


Accumulated air mass distribution for 1 km shells along a 
94° path from 30 km. 


Ray geometry and air mass values for 2 km shells along 91°, 
92° and 93° paths from 30 km. Only one half of the symmet- 
rical air mass distribution is shown. The minimum altitudes 
define the layers for the inversion. 


Integrated absorption distribution of NO for 1 km shells along 
a 93° path from 30 km. 


Integrated absorption distribution of NO 
a 93 path from 30 km. 


2 for 1 km shells along 


Summary of numerical examples for observer altitude of 30 
km (see text for details). 


Air mass distribution for 1 km shells along a 94° path from 
30 km. 
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APPENDIX C 


This appendix includes three sets of tables and figures in this order. 


1) Air mass distributions for balloon altitudes of 15, 20, 25, 30, 35, 
40, 45 and 50 km each for apparent zenith angles of 90. 5°, 91. 0°, 
91.5°, . . . . in steps of 0. 5° until the earth's disk is reached. 


2) Integrated absorption distributions for NO for 30 km and 90. - 
6° WES se 


3) Integrated absorption distributiens for NO, for 30 km and 90. 5°, 
61.0 Ot8, 2 a ES 
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